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ABSTRACT 

A new method for the determination of the extinction in our Galaxy is proposed. 
The method uses surface brightnesses of external galaxies in the B and R bands. 
The observational data have been taken from the ESO-LV galaxy catalogue. As a 
first application of our model we derive the ratio of R band to B band extinction. 
We introduce two methods for computing the ratio which give: 0.62 ± 0.05 (the first 
method) and 0.64 ± 0.06 (the second method) which is in agreement with the recent 
literature value of 0.61. This agreement confirms the validity and efficiency of our 
model and is an independent verification for the standard value of the " total to selective 
extinction" . The method of extinction determination introduced in this paper will be 
explored in subsequent publications. 

Key words: dust, extinction - methods: statistical - Galaxy: general - galaxies: 
fundamental parameters 



1 INTRODUCTION 



Photometric parameters of external galaxies have been used 
several times to derive extinction in our Galaxy. Most au- 
thors use galaxy colours: de Vaucouleurs & Buta (1983), 
Holmberg (1974), Peterson (1970), Sandage (1973), Teeriko- 
rpi (1978). Some of them use galaxy surface brightnesses: 
de Vaucouleurs & Buta (1983), Holmberg (1958) while the 
others use galaxy absolute magnitudes: Peterson (1970), 
Teerikorpi (1978). 

In this paper, we propose a new method of studying 
extinction using surface brightnesses of external galaxies in 
two bands. We apply this method to surface brightnesses in 
B and R bands taken from The Surface Photometry Cata- 
logue of the ESO-Uppsala Galaxies (Lauberts & Valentijn, 
1989, hereafter ESO-LV). The basic idea of the method have 
been shortly described by Choloniewski & Valentijn (1991) 
but the final version of the method described in this paper 
is in some important items different. 

In comparison with earlier methods which use photo- 
metric parameters of external galaxies (see the beginning of 
this Section) our method can be regarded as one which use 
both surface brightnesses of galaxies as well as their colours. 

The expression for the extinction in the B band as used 
in this paper (equation 2) contains a constant, Ar/Ab, the 
ratio of the extinction in the R and the B band. It would 
be easy to take this ratio from the literature and use it for 
extinction determination. But we set out for a more ambi- 
tious task: we try to extract the Ar/Ab ratio from ESO-LV 



surface brightness data alone. An agreement between our 
determination and the literature value would confirm the 
correctness of our model and would demonstrate that the 
surface brightnesses of external galaxies are a good extinc- 
tion indicator. 

We present here two methods for estimating Ar/Ab- 
Both of them rely on the simple notion that extinction 
should not depend on morphological types of the galaxies, 
whose surface brightnesses are used for computing extinc- 
tion. 



2 THE SAMPLE 

The ESO-LV galaxy catalogue contains many different sur- 
face brightness values for every galaxy. Three of them are 
given in two bands (B and R): surface brightness at half to- 
tal B light radius denoted in ESO-LV as )j,f and fi^, average 
central surface brightness within 10 arcsec diameter circular 
aperture denoted as Hq and ^ and central surface bright- 
ness in fit of generalized exponential to octants denoted as 

l l oc 



and ^oct- We use throughout this paper surface bright- 
ness at half total B light denoted hereafter in this paper as 
Hb and fiR. We preferred this measure, as the central sur- 
face brightness might be affected by saturation effects and 
nuclear emission, while the octant surface brightness repre- 
sents extrapolated values. Although the differences of the 
results obtained using these three kinds of surface bright- 
ness values are small, our analysis indicates that the sur- 
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face brightness at half total B light produce results with the 
smallest dispersions. 

We use in this paper the morphological types of galaxies 
listed in the ESO-LV catalogue. However some of them (re- 
ferred in ESO-LV as parametrical ones) are made by means 
of colour indicies which, in turn, are extinction dependent. 
Since it is vital in this paper to use morphological types 
which are absolutely extinction independent we decided to 
reject all galaxies with such morphological types (marked in 
the ESO-LV catalogue with Tfi ag equal to 4). 

Using any astronomical sample for statistical purposes 
one has to specify its definition which clearly describes which 
objects belong to the sample. Generally galaxy samples can 
be apparent diameter-limited (contain galaxies which have 
an apparent diameter larger than a certain diameter limit) 
or magnitude-limited (contain galaxies which are brighter 
than a certain magnitude limit). The ESO-LV catalogue is 
diameter-limited: it contains galaxies which have visual ap- 
parent diameter in the ESO Quick Blue Survey greater or 
equal to 60 arcsec (D org > 60 arcsec). 

However, this completeness limit depends on morpho- 
logical type T. But when we take galaxies which have 
D org > 100 arcsec the catalogue is complete for all morpho- 
logical types (Valentijn 1990, 1994, Huizinga & van Albada 
1992, Impey, Bothun & Malin 1988). So we decided to use 
in this paper galaxies which have D org > 100 arcsec. 

The last (minor) exclusion applied to the original ESO- 
LV sample was the rejection of galaxies which have exces- 
sive (probably wrong) colours: [ib — Hr smaller than and 
greater than 2.1. The exclusion rejects only 18 galaxies. 

The final sample, which has been used for the analysis 
in this paper contains 2783 galaxies, 19 per cent of the full 
original ESO-LV database (which contains 15467 objects). 



3 THE MODEL 

3.1 The formula for extinction 

As one can see in Fig. 1 the surface brightness in the B band 
(/j-b) is strongly linearly dependent on (correlated with) sur- 
face brightness in the R band ([j.r). The amplitude of this 
linear dependence is as high as five magnitudes (similar in 
both bands) and has strictly astrophysical origin (it can not 
be caused by extinction which is for most galaxies substan- 
tially less than one magnitude). 

The Galactic extinction- free values of /is and (j,r 
should, for a given galaxy, increase due to extinction by Ar 
and Ab respectively (the amount of foreground extinction in 
B and R band respectively). This means that on the (/is,/i/j) 
plane extinction "moves" galaxies along the direction with 
slope r: 

Ar m 
r =A~B (1) 

which describes the ratio of extinction in R band to the 
extinction in B band. 

Taking into account the above facts we may expect 
that extinction-free surface brightnesses of galaxies lie on 
a certain straight line. The distance to this zero extinction 
straight line measured in the direction with slope r is pro- 
portional to the Galactic extinction (see Fig. 2) . 

We test our assumption about the linearity of the zero 



extinction line by fitting straight line to the data on the 
(MBjMh) plane (separately for every morphological type T). 
Since the data have errors in both coordinates we have used a 
so called orthogonal fitting procedure (see Feigelson & Babu 
1992). We found no systematic residuals from the resulting 
lines, confirming the validity of our linearity assumption. 

If we describe the slope of the zero extinction straight 
line as s _1 (we use the inverse of s here for symmetry reasons 
- see equations 6, 7, 8, 14 and 15) we can express extinction 
in the B band as: 



(J.B — S flR 



(2) 



1 — r s 

where c is a certain constant. 

Equation (2) gives a ready-to-use formula for extinction 
but in order to use it we have to know three parameters: r, 
s and c. The remaining part of this paper is mainly devoted 
just to obtain these. 

3.2 Selection effects 

The sample used in this paper is apparent diameter (D org ) 
limited. We have checked, that the crucial parameters for 
the present analysis : fiB and fiR are statistically indepen- 
dent of apparent diameter. As we find no dependency this 
implies that selection effects do not influence the values of 
Hb and /j,r, so our formula for Ab, which uses just these two 
observational values is not influenced by selection effects (is 
unbiased). See Choloniewski (1991) for a general discussion 
of photometric selection effects. 

3.3 Relation between r and s 

Let us define now a new variable q as: 

q = fiR- r [i B ■ (3) 

This variable (closely related to the Q parameter defined 
for stars by Sharpless 1963 in a very similar context) is 
extinction-independent because it is equal to the vertical 
distance on the {hb,Hr) plane between a galaxy and extinc- 
tion direction line. As is clearly visible on Fig. 2 the distance 
to the extinction direction line is insensitive to extinction be- 
cause extinction "moves" galaxies on the ([ib,Hr) plane in 
the direction parallel to that line. So, the variable q should 
be statistically independent on As which implies that these 
two variables should be uncorrelated: 

g(A B ,q)=0 (4) 

so their covariance should be equal to zero as well: 

cov(A B ,q) = 0. (5) 

When we replace in equation (5) Ab and q by expressions 
given in equations (2) and (3) respectively and apply the 
law of error propagation (Brandt 1970) we obtain a relation 
which binds the variables r and s: 

ra 2 (p B ) + sa 2 (p R ) = (1 + r s) cov(hb, Hr) (6) 

where a 2 denotes the square of the standard deviation (vari- 
ance). Equation (6) defines on the (hb,(J<r) plane a hyper- 
bola (see Fig. 3) and may be solved with respect to s: 



ra 2 (pL B ) - cov(iib,Hr) 
rcov(p B ,riR) - o- 2 ([ir) 



(7) 
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and r: 



r — 



sct 2 (hr) - covjfiB, (J,r) 



(8) 



3.4 Determination of c 

Our method can provide only relative extinction, namely: 
extinction plus an unknown constant. So, we are free to nor- 
malize our extinction estimate to have: 



< A B >= 0, 



(9) 



what means that our extinction is the relative extinction 
compared to an overall mean. Equation (9) applied to equa- 
tion (2) gives an expression for the constant c: 

< Hb > - s < n R > 



(10) 



3.5 How to compute? 



The average values (more precisely: the expectation values) 
present in equations (9) and (10) can be easily computed 
using the well known formula: 



< x > = 



En 
i = l '■ 



(11) 



where summation is done over the data in the galaxy sample. 
The standard deviations and covariances present in most 
equations in this paper ( can be computed from the following 
formulae: 

a 2 (x) =< x 2 > - <x> 2 (12) 
cov(x,y) —<xy> — <x><y> (13) 
where x and y denotes any variable used in this paper. 

3.6 Statistical and geometrical context 

Equations (2) and (3) define in fact a linear transformation 
of the "old" variables \ib and \h,r into the "new" variables 
Ab and q. The transformation is defined in such a way that 
the "new" variables have to be uncorrelated. The coefficients 
of such a linear transformation are bound by equation (6). 

We solve here a general problem: how to transform lin- 
ear "old" variables x and y into the "new" ones X and Y: 



X — x — sy 
Y = y — rx 

while X and Y are uncorrelated: 
g(X,Y) = 



(14) 
(15) 

(16) 



The solution of this problem, lies in the condition for s and 
r: 



ro~ 2 {X) + so~ 2 {Y) = (1 +rs)cov(X,Y) 



(17) 



Equation (6) is of course the special case of this equation. 

Geometrically, the linear transformation of two "old" 
variables into the "new" ones is equivalent to applying a 
certain new coordinate system, which uses two "new" axes: 
the X "new" axis is tilted to the "old" x axis by a slope 
r, while the "new" Y axis is tilted to the "old" y axis by a 
slope s. 



Our solution is closely related to the procedure which 
leads to the so called orthogonal regression line (see Feigel- 
son & Babu 1992). Such a line is constructed by rotating 
coordinates in a way to obtain uncorrelated data in "new" 
coordinates. Our solution can be regarded as a generaliza- 
tion of this procedure: we can easily reconstruct an orthogo- 
nal regression line by substituting in equation (17) s = —r. 



4 DETERMINATION OF A R /A B : FIRST 
METHOD 

We compute in this Section the ratio of extinction in R band 
to the extinction in B band denoted as r (see equation 1) 
using the simple fact that extinction As must be statistically 
independent on galaxy morphological type T. 



4.1 Mathematical approach 

The statistical independence of Ab and T imply that: 
g(A B ,T) = (18) 



which gives: 
cov(A B ,T) = 0. 



(19) 



Taking into account equation (2), using the law of error 
propagation (Brandt 1970) and assuming that the parame- 
ters s and c do not depend on T we have: 



cov(p B ,T) 
cov(p,R,T) ' 



(20) 



This value for s substituted into equation (8) gives the for- 
mula for r: 



cov(hb,T) 2, s 
cov(n R ,T) " \PR) 



cov(hb,^r) 



Z(mh,'t) cov(^b,Vr.) - o- 2 (p B ) 



(21) 



Equations (20) and (21) hold under the assumption that 
s and c do not depend on T. Let us check now whether 
this is true by computing s and c, for different values of r 
using equations (7) and (10). The uncertainties for s and c 
have been computed using a so called jacknife method (see 
Quenouille 1956, Tukey 1958 and Efron & Tibshirani 1986). 
The results are given on Fig. 4 and 5. It is evident that s and 
c are T dependent, but for morphological type T between 
2.5 and 6.5, both parameters are constant for any value of 
r in the limits of accuracy achieved here. 

For this reason we computed r using equation (21) tak- 
ing only those galaxies for with T is greater than 2.5 and 
less than 6.5. The resultant subsample contains 1290 objects 
(approximately half of the whole sample). We computed for 
these galaxies all the components of the equation (21): 



cov(lib,T) = 279.032 
cov(hr,T) = 416.098 
a 2 (fi B ) = 0.432634 
CT 2 (/i fl ) = 0.530865 
cov(p B ,tiR = 0.440979 
What finally gives: 
r = 0.62 ±0.05 



(22) 
(23) 
(24) 
(25) 
(26) 

(27) 
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where the uncertainty of the result was again computed us- 
ing the jacknife method. 

4.2 Graphical approach 

The idea which has been used here for computing r (A B 
does not depend on T) can be applied in another way. Let 
us compute the average extinction (according to equations 
2 and 7) and plot it as a function of morphological type 
(denoted as < Ab >t) for different values of r - see Fig. 6. 
It is clearly visible that only for r — 0.6 the extinction Ab is 
not a function of morphological type T. So we arrive again 
at a similar result as expressed in equation (21). 



5 DETERMINATION OF A R /A B : SECOND 
METHOD 

We compute in this Section the ratio of extinction in R band 
to the extinction in B band denoted as r (see equation 1) 
using the fact that the standard deviation of extinction Ab 
must be the same for every galaxy morphological type T. 

According to equation (2) and applying the law of er- 
ror propagation (Brandt 1970) the value of a 2 (As) can be 
expressed as: 

2, A \ ct 2 (mb) + s 2 a(p R ) -2scov{ij,b,Hr) . os n 
° {Ab) = (1 -rsY (28) 

When we substitute s present in equation (28) by expres- 
sion given in equation (7) the resultant equation for o 2 (Ab) 
depends only on r. When we use it for two different mor- 
phological types (say Ti and Tb) we have: 



CTTi {A t 



(TT 2 {Al 



(29) 



which can be solved with respect to r. Such solution is pre- 
sented graphically on Fig. 7 for Ti = —3 and T-2 = 5. The 
appropriate two curves intersect at r = 0.66 which is the 
solution for r for this particular pair of T. 

The values of ot{Ab)t as a function of r for all morpho- 
logical types (T = -5, ... , 10) are presented in Fig. 8. The 
curves intersect in many different points so we have many 
different "solutions" for r. This motivated us to compute the 
relative standard deviation of <7T(AB)r, namely: 



D(r) 



cr(aT{A B )r) 

< T {A B )r > 



(30) 



The minimum of the function D(r) will be the solution for 
r. For two morphological types: T=-5 and T=9 aT(As) r 
strongly differs from the mean so we decided to omit these 
two types in our final computations in this Section (the re- 
jected curves are marked on Fig. 8 by a dotted lines). 

The reason for existence of such two "outliers" is prob- 
ably that thay have different intrinsic extinction-free scat- 
ter with respect to the "zero extinction" straight line on the 
(ij*b,IJ.r) plane than for other galaxies: this scatter is smaller 
(in comparison with most galaxies) for T=-5 and greater for 
T=9. 

The function D(r) has a minimum for: 



0.64 ±0.06 



(31) 



where the uncertainty has been computed by the jacknife 
method (with respect to T). 



6 DISCUSSION 

The second method for determination of the An/ As ratio 
works only due to the fact that the parameter s is T depen- 
dent. Moreover, the stronger this dependence is, the more 
accurate are the results given by the second method. The 
opposite situation is for the first method: it works only for 
such galaxies for which the s parameter is not T dependent. 
So the first method can survive only due to galaxies with 
the same s parameter which is in contradiction to the sec- 
ond method which feeds itself by the differences in s. While 
complementary in nature, both tests gives similar results. 

Both methods have their specific disadvantages. The 
first one does not use the whole sample. The second one 
relies on the assumption that the extinction-free scatter of 
galaxy surface brightnesses is the same for every morpho- 
logical type (what is not always true). 

The value for Ar/Ab which was recently published is 
0.61 (Schlegel, Finkbeiner & Davis 1998). It takes into ac- 
count several observational and instrumental factors and is 
based on formulae published by Cardelli, Clayton & Mathis 
1989 and O'Donnell 1994. Our results: 0.62 ± 0.05 (the first 
method) and 0.64 ± 0.06 (the second method) are in very 
good agreement with this value. It confirms the correctness 
of the model introduced in this paper and practically demon- 
strates that two band surface brightness data of external 
galaxies are a good extinction indicator. 

Since Ar/Ab ratio is closely connected to "total to se- 
lective extinction ratio" Av/EiB — V) our result confirms 
its standard value which is equal to 3.1. 

The method applied in this paper to obtain Ar/Ab ra- 
tio relies on the assumption that extinction does not depend 
on morphological type. In future applications to other data 
(e.g. SDSS) one can use, instead of morphological type, other 
extinction independent quantities like: axial ratio, redshift 
or effective diameter. 
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Figure 1. Surface brightnesses at half total B light radius of galaxies in the B and R bands for four different morphological type T 
intervals. The data are from a subsamplc of the ESO/LV galaxy catalogue which has been used for all the computations presented in 
this paper. 
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Figure 2. The idea of extinction measurements using surface brightness of galaxies in two bands: fig and fiR. The distance of a galaxy 
(marked as an asterix) to the zero extinction line measured parallel to the extinction direction line is proportional to the Galactic 
extinction. The q parameter is extinction independent. The slope of the zero extinction line is s _1 while the slope of the extinction 
direction line is r = An/Ag. 
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Figure 3. The relation between the parameters s and r = Ar/Ab derived from the requirement that extinction Ab and the parameter 
q are statistically independent (as described in equation 6). 




Figure 4. The parameter s as a function of morphological type T for various values of r = An/Ag. Error bars represent standard 
deviation (la). 
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Figure 5. The parameter c as a function of morphological type T for various values of r = Ar/Ab- Error bars represent standard 
deviation (la). 
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Figure 6. The average extinction in B band < Ag > as a function of morphological type T for various values of r = An/Ag. The 
extinction should not depend on morphological type which is true only for r=0.6. Error bars represent standard deviation (la). 
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Figure 7. The standard deviation of extinction in B band <t(Ab) as a function of r = Ab/Ab for galaxies with T = —3 and T = 5. The 
point of intersection of these two curves at r=0.66 gives the solution for Ar/Ab since cf(Ab) should be the same for any T. 
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Figure 8. The same as Fig. 7 but for all morphological types T= -5, ... , 10). The dotted line refers to galaxies with T=-5 (at the 
bottom) and T=9 (at the top). These galaxies strongly differ from the mean so they have been rejected in final computations. 
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Figure 9. The function D(r) which describes the relative standard deviation of the <j(Ab)- The minimum of this function for r=0.64 
± 0.06 represents the solution for r = Ar/Ab since a(Ag) should be as close as possible to each other (in ideal case the same) for all 
morphological types. 



